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We apply standard post-Newtonian methods in general relativity to locate the 
innermost circular orbit (ICO) of irrotational and corotational binary black-hole 
systems. We find that the post-Newtonian series converges well when the two 
masses are comparable. We argue that the result for the ICO which is predicted 
by the third post-Newtonian (3PN) approximation is likely to be very close to the 
"exact" solution, within 1% of fractional accuracy or better. The 3PN result is 
also in remarkable agreement with a numerical calculation of the ICO in the case 
of two corotating black holes moving on exactly circular orbits. The behaviour of 
the post-Newtonian series suggests that the gravitational dynamics of two bodies 
of comparable masses does not resemble that of a test particle on a Schwarzschild 
background. This leads us to question the validity of some post-Newtonian resum- 
mation techniques that are based on the idea that the field generated by two black 
holes is a deformation of the Schwarzschild space-time. 

1 Introduction 

The "standard" post-Newtonian approximation, or expansion when the speed 
of light c — > +oo, is at the basis of an important body of research, which has 
provided us in the past with our best picture of the motion of compact objects 
in general relativity. We quote the pionneering works of EinsteirP, Droste^, 
and DeSitteJ^l the landmark analysis due to Einstein, Infeld and Hoffmanr0of 
the dynamics of N separated bodies at the first post-Newtonian orde r (1PN , or 
1/c 2 ), and the seminal papers by Chandrasekhar and collaboratora™^ con- 
cerning the equations of motion of extended fluid systems, up to the 2.5PN 
level (at which order appears the first effect due to the reaction to the emis- 
sion of gravitational radiation). In the case of two compact objects (neutron 
stars or blac k holes ), we possess the 2.5PN equations of motion of the bi- 
nary pulsar^ jj-j 1 y 1 1 ~H and the 3P N eq uat ions of motion of inspiralling compact 
binarie d ' ' ' ' ' < | io | ib> | zu | /i | Regarding the problem of the gravita- 
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tional radiation emitted by inspiralling compact binaries, we hav e under con- 
trol most of the gravitational-wave flux up to the 3.5PN or 
including the specific effects of wave taila^EZI 

In this paper we focus our attention on the question of the dynamics of 
black-hole binary systems (henceforth we assume that the compact objects 
are black holes) and its recent resolution up to the 3PN approximation, corre- 
sponding to the formal order 1 /c 6 beyond the Newtonian force law. On the one 

hand the ADM-Hamiltonian formalism of general relativity has been applied 

11 011 L \\ 

at the 3PN order by Jaranowski and Schafe r ' ° \ and Damour, Jaranowski 
and SchafeOT On the other hand a direct 3PN iteration of the equations 
of motion — instead of a Hami lton ian — in harmonic coordinates (extend- 
ing the method proposed in has been implemented by Blanchet and 
and de Andrade, Blanchet and These two inde- 
pendent approaches have succeeded; it has been shown that there exists a 
unique transformation of the particle's dynamical variables that changes the 
3PN harmonic-coordinates Lagrangian^ 1 into a Lagrangian whose Legendre 
transform is exactly identical to the 3PN ADM-coordinates Hamiltonian^^. 

In the previous approaches the two black holes are modelled by point 
particles, solely described by their masses m\ and mi. This is consistent with 
the very spirit of the post-Newtonian method. The point-particle description 
has to be supplemented by a process of regularization of the self field of point 
parti cles. The standard r egularization a la Hadamard is at the basis of the 
workJ^ -^-^l ^'fl^W ^ I (an extended version of this regularization has been 
defined in Refs™). Unfortunately it was shown that at the 3PN order 
the equations of motion of black holes contain a particular coefficient, i.e. 
the "static" ambiguity uj s in the ADM-Hamiltonian formalism^ W W or 
the parameter A in the harmonic-coordinates approao < l 1 g l i v l z v l /1 P i which 
cannot be fixed by the Hadamard regularization. The value of this coefficient 
has been obtained by means of a dimensional regulari zati on instead of the 
Hadamard one within the ADM-Hamiltonian formalism^!. We shall discuss 
some implications of this result (i.e. uj s = 0) regarding the validity of the 
post-Newtonian expansion. 

Let us look at the conserved energy of the black-hole binary in the center- 
of-mass frame at the 3PN order. Technically this energy follows from the 
3PN harmonic-coordinates Lagrangian or equivalently from the 3PN ADM- 
coordinates HamiltoniarPSI. The energy is conserved only when we neglect the 
radiation reaction damping at the 2.5PN order. The time derivative of the 



a See Eq. JoJ below for the relation linking together A and uj s 
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energy is equal to the radiation reaction effect, but for the present discussion 
we are interested only in the conservative part of the dynamics which is com- 
posed of the Newtonian, 1PN, 2PN and 3PN approximations. Specializing 
to the case of orbits which are circular (apart from the gradual radiation- 
reaction inspiral) 6 yields then the center-of-mass energy E at the 3PN order 
for circular orbits as a function of the frequency u> of the orbital motion. 

Here we want to assess the validity of the post-Newtonian approxima- 
tion. More precisely we address, and to some extent we answer, the following 
questions. How accurate is the post-Newtonian expansion for describing the 
dynamics of binary black hole systems ? Is the innermost circular orbit (ICO) 
of binary black holes, defined by the minimum of the energy function E(u), 
accurately determined at the highest currently known post-Newtonian order ? 
This question is pertinent because the ICO represents a point in the late stage 
of evolution of the binary which is very relativistic (orbital velocities of the 
order of 50% of the speed of light). How well does the 3PN approximation as 
compared with the prediction provided by numerical relativity ? What is the 
validity of the various post-Newtonian resummation techniqueJ^SEHl which 
aim at "boosting" the convergence of the standard post-Newtonian approxi- 
mation ? 

The previous questions are very interesting but difficult to settle down rig- 
orously. Indeed the very essence of an approximation method is to cope with 
our ignorance of the higher-order terms in some expansion, but the higher- 
order terms are precisely the ones which would be needed for a satisfying an- 
swer to these problems. So we shall be able to give only some educated guesses 
and/or plausible answers, that we cannot justify rigorously, but which seem 
very likely from the standard point of view on the post-Newtonian theory, in 
particular that the successive orders of approximation get smaller and smaller 
as they should (in average), with only few accidents occuring at high orders 
where a particular approximation would be abnormally large with respect to 
the lower-order ones. Admittedly, in addition, our faith in the estimation we 
shall give regarding the accuracy of the 3PN order for instance, comes from 
the historical perspective, thanks to the many successes achieved in the past 
by the post-Newtonian approximation when confronting the theory and ob- 
servations (e.g. of the Solar system dynamics and the binary pulsar). It is 
indeed beyond question, from our past experience, that the post-Newtonian 
method does work. 

There are many other related issues that we shall not address. For in- 



b We know that most inspiralling compact binaries have circular orbits because of radiation- 
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stance : is the best presently known post-Newtonian approximation, i.e. 3PN, 
sufficient for the problem of the coalescence of two black holes in a foresee- 
able future c ? A related problem is the qualitative and quantitative influence 
of the high-order radiation reaction effects which have to be superposed to 
the conservative part of the dynamics. Finally we shall discuss only the bi- 
nary's equations of motion, and leave aside the problem of the radiation field. 
Certainly the accuracy of the 3.5PN gravi tationa l-radiation flux of black-hole 
binaries which has been computed in Refs.EMSl should be discussed in a sim- 
ilar way. 

Basically, the point we would like to emphasize 11 is that the post- 
Newtonian approximation, in stand ard form (without using the resummation 
techniques advocated in RefsPSESI) ; i s a bi e to located the ICO of two black 
holes with an excellent accuracy of the order of 1%, and perhaps much better 
than that. At first sight this statement is rather surprising, because the dy- 
namics of two black holes at the point of the ICO is so relativistic. Indeed one 
often ears about the "bad convergence" , or the "fundamental breakdown" , of 
the post-Newtonian series in the regime of the ICO. However our estimates 
do show that the 3PN approximation is very good in this regime, and they 
are also confirmed by the remarkable agreement, we shall detail below, with 
a numerical calculation by Gourgoulhon et aU ' 'oi the ICO in the case of 
two black holes moving on exactly circular orbits ("helical symmetry"). When 
comparing the post-Newtonian approximation with the numerical simulation 
we face an interesting problem : since the numerical work has been done for 
corotating black holes, which spin with the orbital velocity ui, the effects of 
spins, appropriate to two Kerr black holes, are to be taken into account within 
our post-Newtonian framework. 

Another point we shall develop is that most probably the general- 
relativistic dynamics of two objects with comparable masses is qualitatively 
different (in addition of being far more complicated), than the dynamics of 
a "test" particle with geodesic motion on a fixed background Schwarzschild 
space-time. Our argument has something to do with the value lo s — ob- 
tained in Ref for the 3PN regularization ambiguity parameter. Indeed this 
value, if we believe that it is really the prediction of general relativity, sug- 
gests that the two-body interaction is not "Schwarzschild-like" , in the sense 
that it does not seem to admit a light-ring singularity similar to the one of 



c We have in mind the 3PN approximation considered as an input for the numerical calcu- 
lation of the black-hole coalescence (assuming that a clear method would exist for imple- 
menting the post-Newtonian initial conditions into a numerical scheme). 
We are following the recent study in 
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the Schwarzschild space-time. We then argue that this fact sheds a doubt 
on the validity, at high post-Newtonian orders, of some resummation tech- 
niques (like Pade approximanta2SI) that are based on the assumption that the 
field generated by two bodies of comparable sizes is a "deformation" of the 
Schwarzschild metric. A contrario we shall find that the value o; s = gives 
some convincing evidence that the standard post-Newtonian approximation, 
based on Taylor rather than Pade approximants, is very accurate. 

2 The binding energy of two black-holes 

The binding energy, in the center-of-mass frame, is defined as the invariant 
energy associated with the conservative part of the binary's 3PN dynamics 
(we ignore the radiation reaction effect at the 2.5PN order). Restricting our 
consideration to circular orbits, the energy is a function of a single variable, 
which can be chosen to be the distance r between the two particles in a given 
coordinate system, or, better, the frequency u> — of the orbital motion (P 
is the orbital period) . We introduce for convenience the particular frequency- 
related parameter 



The mass-energies of the black holes are mi and m 2 (they take notably into 
account the rotational energies); the total mass is M = mi + to 2 . It is 
important to express the energy in terms of the frequency-related parameter 
x, instead of some coordinate distance r, because the function E(x) then takes 
an invariant expression (the same in different coordinate systems). 

Having in hands the circular-orbit energy, we then define the innermost 
circular orbit (ICO) as the minimum, when it exists, of the energy function 
E(x). T he defin ition is motivated by our comparison with the numerical 
calculatio n^ ! ^ ! because this is precisely that minimum which is computed 
numerically. In particular, we do not define the ICO as a point of dynamical 
general-relativistic unstability. The circular-orbit energy, developed to the 
3PN order, is of the form 



The first term is the rest-mass, the next one, proportional to x, is 
the Newtonian term, and then we have many post-Newtonian correc- 




(1) 
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tions, the coefficients of which are known at present only up to the 3PN 

0Td Jm-mmmwm, and given by 

ai(u) =--- — , (3) 
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We make use of the useful ratio between the reduced and total masses 



fi mim 2 
v = — where u = — — — . (o) 
M M w 

This ratio is interesting because of its range of variation, 

< v < \ , (7) 

where v = j in the equal-mass case and v — > in the test-mass limit for one 
of the bodies. We shall investigate the value of the ICO as predicted by the 
3PN energy J5J- 101 in Section 4. 

The 3PN coefficient a^{y) i nvolves the regularization-ambiguity param- 
eter A introduced in Ref s! ' ' and due to a physical incompleteness in the 
Hadamard methoc P^^ 9 ! for regularizing the self-field of poi nt particles. This 
p arameter is equivalent to the parameter to s in Refsl ^ l ^ and related to it 

x 3 1987 (Z\ 

A = -TT w - " 3080- (8) 

It has been argued in that the numerical value of w s could be ~ —9, 

because for such a value some different rcsummation techniques, when they 
are implemented at the 3PN order , gi ve approximately the same result for the 
ICO. Even more, it was suggested^ that oj s might be precisely equal to w*, 
with 
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-9.34- •• . (9) 
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However, a more recent computation performed with the help of a dimensional 
regularization instead of the Hadamard regularization, within the ADM- 
Hamiltonian formalism^! has yielded 



— — < io > 

Here we adopt lu s = as the "correct" value predicted by general relativity for 
the ambiguity parameter. Note that the result (|1(JH is quite different from uj* 
given by Eq. : this suggests, accor ding to Ref! 4 ™ that different resumma- 
tion techniques, viz Pade approximant inland effective-one-body methodJ23 
which are designed to "accelerate" the convergence of the post- Newtonian se- 
ries, do not in fact converge toward the same exact solution (or, at least, not 
as fast as expected). 

The appearance of one ambiguity parameter at the 3PN order is inter- 
esting in con nection with the so-called effacing principle satisfied by general 
relativity^, according to which the equations of motion and the radiation field 
of gravitationally condensed objects should depend only on their relativistic 
masses and not on the detailed features of their internal structure — we ne- 
glect the tidal effects between the objects. In general relativity this principle 
follows from the strong equivalence principle (which differs from the Einstein 
equivalence principle by the inclusion of bodies with self-gravitational interac- 
tions and of experiments involving gravitational forces). Indeed, in the freely 
falling frame of one of the (spherically symmetric) compact objects, we can 
apply the Birkhoff theorem and find that the external vacuum field depends 
only on the mass. In consequence we should expect that the parameter A is a 
pure number (e.g. a rational fraction), independent of the internal structure 
of the compact objects. 

It would be interesting to confirm Eq. (|f 0|) by an independent calculation, 
hopefully without resort to any regularization scheme. We have in mind a cal- 
culation valid for extended ("fluid") systems, taking a priori into account the 
internal structure of the fluids. Such a calculation would provide, in addition 
to the determination of A, an explicit verification of the effacing principle of 
general relativity. Notice that from the point of view of a calculation valid 
for extended fluids, it is hard to believe that the ambiguity parameter could 
depend on tt 2 like in Eq. ©• 

In the limiting case v — > 0, the energy 0-© reduces to the 3PN ap- 
proximation of the exact expression for a test particle in the Schwarzschild 
background, 
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(11) 



The minimum of that function or ICO occurs at x^qq = ^, and we have 
£f£h = ^c 2 - 1) . We know that the ICO in the case of the Schwarzschild 

metric is also an innermost stable circular orbit (or ISCO), i.e. it corresponds 
to a point of dynamical unstability. Another important feature of Eq. is 
the singularity at the value ^fight-ring = § which corresponds to the famous 
circular orbit of photons in the Schwarzschild metric ( "light-ring" singularity) . 
This orbit can also be viewed as the last unstable circular orbit. We can check 
that the post-Newtonian coefficients a^ ch = a n (0) corresponding to Eq. 
are given by 

Sch 3"(2n-l)!!(2n-l) 
G ™ ~ 2»(n + l)l ' [ ' 

They increase with n by roughly a factor 3 at each order. This is simply the 
consequence of the fact that the radius of convergence of the post-Newtonian 
series is given by the Schwarzschild light-ring singularity at the value |. We 
may therefore recover the light-ring orbit by investigating the limit 

Sch i 

lim n ~ 1 - - - r Sch fnl 

„^+oo a Sch _ 3 _ ^light-ring ■ K 1 ^) 



3 Accuracy of the post-Newtonian expansion 

Let us now discuss the accuracy of the 3PN approximation when estimating 
the ICO of binary black holes. First of all we make a few order-of-magnitude 
estimates. At the location of the ICO we shall find (see the next Section) that 
the frequency-related parameter x defined by Eq. (:1) is approximately of the 
order of 20%. Therefore, we might a priori expect that the contribution of the 
1PN approximation to the energy at the point of the ICO should be of that 
order. For the present discussion we take the pessimistic view that the order of 
magnitude of an approximation represents also the order of magnitude of the 
higher-order terms which are neglected. We see that the 1PN approximation 
should yield a rather poor estimate of the "exact" result, but this is quite 
normal at this very relativistic point where the orbital velocity is - ~ x 1 / 2 ~ 
50%. By the same argument we infer that the 2PN approximation should do 
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much better, with fractional errors of the order of x 2 ~ 5%, while 3PN will 
be even better, with the accuracy x 3 ~ 1%. 

The simple order-of-magnitude estimate suggests therefore that the 3PN 
order should be close to the "exact" solution for the ICO to within 1% of 
fractional accuracy. We think that this is very good, and we should even 
remember that this estimate is pessimistic, because we can reasonably expect 
that the neglected higher-order approximations, 4PN and so on, are in fact 
smaller numerically (e.g. of the order of x 4 ~ 0.2%). But let us keep for the 
present discussion the 1% guess for the accuracy of the 3PN approximation. 

Now the previous estimate makes sense only if the numerical values of 
the post-Newtonian coefficients in Eqs. ©-© stay roughly of the order of 
one. If this is not the case, and if the coefficients increase dangerously with 
the post-Newtonian order n, one sees that the post-Newtonian approximation 
might in fact be very bad. It has often been emphasized in the litterature 
(see e.g. Refs PESEHI ) that in the test-mass limit v — ► the post-Newtonian 
series converges slowly, so the post- Newtonian approximation is not very good 
in the regime of the ICO. Indeed we have seen that when v — the radius of 
convergence of the series is | (not so far from x\^ — |), and that accordingly 
the post-Newtonian coefficients increase by a factor ~ 3 at each order. So it 
is perfectly correct that in the case of test particles in the Schwarzschild 
background the post-Newtonian approximation is to be carried out to a high 
order in order to locate the turning point of the ICO. 

Let us immediately remark that this negative conclusion does not matter : 
indeed we shall never use the post-Newtonian approximation when v — > 
simply because we know the exact result which is given by Eq. (8) e . Therefore 
we should not worry about the poor convergence of the post-Newtonian series 
in the test-mass limit. The post-Newtonian method is useless and even one 
might say irrelevant when considering the motion of a test particle around a 
Schwarzschild black hole. 

What happens when the two masses are comparable {y = 4) ? It is clear 
that the accuracy of the post-Newtonian approximation depends crucially on 
how rapidly the post-Newtonian coefficients increase with n. We have seen 
that in the case of the Schwarzschild metric the latter increase is in turn related 
to the existence of a light-ring orbit. For continuing the discussion we shall 
say that the relativistic interaction between two bodies of comparable masses 
is "Schwarzschild-like" if the post-Newtonian coefficients a n (j) increase when 
n — > +oo. If this is the case this signals the existence of something like a light- 
ring singularity which could be interpreted as the deformation, when the mass 



e The exact result for the radiation field is also known, albeit only numerically. 
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Table 1. Numerical values of the sequence of coefficients of the post-Newtonian series com- 
posing the energy function I2*t-l5t. 





Newtonian 




a.2 (y) 




v = 


1 


-0.75 


-3.37 


-10.55 


" = i < - 


-9.34 1 


-0.77 


-2.78 


-8.75 



v=j lo s = 1 -0.77 -2.78 -0.97 



ratio v is "turned on", of the Schwarzschild light-ring orbit. By analogy with 
Eq. (jT3|) we can estimate the location of this "pseudo-light-ring" orbit by 

~ flight-ring^) with n = 3 . (14) 

Here n = 3 is the highest known post-Newtonian order. If the two-body 
problem is "Schwarzschild-like" then the right-hand-side of Eq. I|14|l is small 
(say around g), the post-Newtonian coefficients typically increase with n, and 
most likely it should be difficult to get a reliable estimate by post-Newtonian 
methods of the location of the ICO. So we ask : is the gravitational interaction 
between two comparable masses Schwarzschild-like ? 

In TableHwe present the values of the coefficients a n (v) in the test-mass 
limit v = [see Eq. 1)12(1 for the analytic expression], and in the equal-mass 
case v = j when the ambiguity parameter takes the "uncorrect" value oj* 
[defined by Eq. @] and the correct one u> s = predicted by general relativity. 
When v = we clearly see the expected increase of the coefficients by roughly 
a factor 3 at each step. Now when u = -r and ui s = uj* we notice that the 
coefficients increase approximately in the same manner as in the test-mass 
case v — 0. This indicates that the gravitational interaction in the case of 
u* looks like that in a one-body problem. The associated pseudo-light-ring 
singularity is estimated using Eq. ljl"H) as 

Zlight-ring(3,W*) ~ 0.32 . (15) 

The pseudo-light-ring orbit seems to be a very small deformation of the 
Schwarzschild light-ring orbit given by Eq. (|13f) . In this Schwarzschild-like 
situation, we should not expect the post-Newtonian series to be very accurate. 

Now in the case v = \ but when the ambiguity parameter takes the 
correct value lo s = 0, we see that the 3PN coefficient 03(2) is of the order 
of minus one instead of being ~ —10. This suggests (unless 3PN happens to 
be quite accidental) that the post-Newtonian coefficients in general relativity 
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do not increase very much with n. This is an interesting finding because 
it indicates that the actual two-body interaction in general relativity is not 
Schwarzschild-like. There does not seem to exist something like a light-ring 
orbit which would be a deformation of the Schwarzschild one. Applying Eq. 
(|14l) we obtain as an estimate of the "light-ring" : 

zii g ht-ring(3,G.R.)~2.88. (16) 

It is clear that if we believe the correctness of this estimate we must conclude 
that there is in fact no notion of a light-ring orbit in the real two-body problem. 
Or, one might say (pictorially speaking) that the light-ring orbit gets hidden 
inside the horizon of the final black-hole formed by coalescence. Furthermore, 
if we apply Eq. (|14J) using the 2PN approximation n ~ 2 instead of the 
3PN one n = 3, we get the value ~ 0.28 instead of Eq. Q16p. So at the 
2PN order the metric seems to admit a light ring, while at the 3PN order it 
apparently does not admit any. This erratic behaviour reinforces our idea that 
it is meaningless (with our present 3PN-based knowledge, and untill fuller 
information is available) to assume the existence of a light-ring singularity 
when the masses are equal. 

It is impossible of course to be thoroughly confident about the validity 
of the previous statement because we know only the coefficients up to 3PN 
order. Any tentative conclusion based on 3PN can be "falsified" when we 
obtain the next 4PN order. Nevertheless, we feel that the mere fact that 
a 3(|) = — 0.97 in Table is sufficient to motivate our (tentative) conclusion 
that the gravitational field generated by two bodies is more complicated than 
the Schwarzschild space-time. This appraisal should look cogent to relativists 
and is in accordance with the author's respectfulness of the complexity of the 
Einstein field equations. 

We want next to comment on a possible implication of our conclusion 
as regards t he so-c alled post-Newtonian resummation techniques, i.e. Pade 
approximant J2BES1 which aim at "accelerating" the convergence of the post- 
Newton ian ser ies in the pre-coalescence stage, and effective-one-body (EOB) 
methoda^2E21 which attempt at describing the late stage of the coalescence of 
two black holes. These techniques are based on the idea that the gravitational 
two-body interaction is a "deformation" — with v < j being the deformation 
parameter — of the Schwarzschild space-time. The Pade approximants are 
valuable tools for giving accurate representations of functions having some 
singularities. In the problem at hands they would be justified if the "exact" 
expression of the energy [whose 3PN expansion is given by Eqs. ©-©] 
would admit a singularity at some reasonable value of x (e.g. < 0.5). In 
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the Schwarzschild case, for which Eq. I|13|) holds, the Pade series converges 
rapidlj®' : the Pade constructed only from the 2PN approximation of the 
energy — keeping only af ch and af ch — already coincide with the exact result 
given by Eq. On the other hand, the EOB method maps the post- 

Newtonian two-body dynamics (at the 2PN or 3PN orders) on the geodesic 
motion on some effective metric which happens to be a j/-deformation of the 
Schwarzschild space-time. In the EOB method the effective metric looks like 
Schwarzschild by definition, and we might expect the two-body interaction to 
own some Schwarzschild-like features. 

Our comment is that the validity of these post-Newtonian resummation 
techniques is questionable because the value lu s = suggests that the two-body 
interaction is not Schwarzschild-like. Our doubt is confirmed by the finding of 
RefP^ (already alluded to above) that in the case of the "wrong" ambiguity 
parameter oj* ~ —9.34 the Pade approximants and the EOB method at the 
3PN order give the same result for the ICO. From the previous discussion we 
see that this agreement is to be expected because a deformed light-ring sin- 
gularity seems to exist with lu*. By contrast, in the case of general relativity 
(lu s = 0), the Pade and EOB methods give quite different results (cf. the 
figure 2 in Ref.^). Such a disagreement, we argue, is due to the fact that the 
assumptions underlying the various resummation techniques are probably not 
fulfilled, so they may converge toward different solutions. Another confirma- 
tion comes from the light-ring singularity which is det erm ined from the Pade 
approximants at the 2PN order [see Eq. (3.22) in as 



£iight-rmg(|, Pade) ~ 0.44. (17) 

This value is rather close to Eq. i|15[l but strongly disagrees with Eq. i|16|) . 
Our explanation is that the Pade series converges toward a theory having 
uj s ~ lo* and so which is different from general relativity. 

Finally we come to the good news that, if really the (absolute value of the) 
post-Newtonian coefficients when v = j stay of the order of one as it seems to, 
this means that the standard post-Newtonian approach, based on the standard 
Taylor approximants, is probably very accurate. The post-Newtonian series 
is likely to "converge well" , with a "convergence radius" of the order of one-^ . 
Therefore the order-of-magnitude estimates we proposed at the beginning of 

* Actually, the post-Newtonian series could be only asymptotic (hence divergent), but nev- 
ertheless it should give excellent results provided that the series is truncated near some 
optimal order of approximation. In this discussion we assume that the 3PN order is not 
too far from that optimum. 
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this Section are probably correct. In particular the 3PN order should be close 
to the "exact" solution even in the regime of the ICO. 

Compare the situation with the case of Schwarzschild, for which we have 
seen that xf^Q — | is rather close to the value of the convergence radius of 
the series given by ^fight-ring = I' nence the poor convergence of the post- 
Newtonian expansion at the location of the ICO. By contrast, when the two 
masses have the same size, we have xico(;j) ~ 0.2 (see the Figure^) which is 
quite far from the "convergence radius" of the series, that we argued is likely 
to be of the order of one. We thus expect the post-Newtonian expansion (in 
the case of two black holes) to be well appropriate near the ICO. 



4 Comparison with a numerical simulation 

We confront the prediction of the standard (Taylor-based) post-Newtonian 
approach with a recent re sult o f numerical relativity by Gourgoulhon, 
Grandclement and BonazzolcP^EH. These authors computed numerically the 
energy of binary black holes under the assumptions of conformal flatness for 
the spatial metric and of exactly circular orbits. The latter restriction is im- 
plemented by requiring the existence of an "helical" Killing vector, time-like 
inside the light cylinder associate d with the circular motion and space-like 
outside. In the numerical approacn '"' L ' v^l there are no gravitational waves, the 
field is periodic in time, and the gravitational potentials tend to zero at spa- 
tial infinity within a restricted model equivalent to solving five out of the ten 
Einstein field equations. Consid ering a n evolutionary sequence of equilibrium 
configurations Gourgoulhon et a/! ' ' obtained numerically the circular-orbit 
energy E(u>) and looked for t he inne rmost circular orbit or ICO. 

The numerical calculation 1 ^^ ! has been performed in the case of coro- 
tating black holes, which are spinning with the orbital angular velocity ui. For 
the comparison we must therefore include within our post-Newtonian formal- 
ism the effects of spins appropriate to two Kerr black holes rotating at the 
orbital rate u>. The importance of the effect of spins in corotating systems of 
neutron stars, for which the ICO is usually determined by the hydrodynamical 
instability rather than by the effect of general relativity, is well knowrPSI We 
expect that these effects play some role in the case of black holes as well. 

The total relativistic mass of the Kerr black hole is given by 9 



S 2 

m 2 = m 2 m + ^-g- , (18) 



9 In the formulas fTgl-iSSt we pose G = 1 
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irr 



where S is the spin, related to the usual Kerr parameter by S = ma, and m 
is the irreducible mass given by mj rr = ^0- (A is the hole's surface area). The 
angular velocity of the corotating black hole is uj = ^S- hence, from Eq. 118|) . 



OS 

S 



us 



2m 3 


1 + 


\/l - ^ 






V m 



(19) 



Physically this angular velocity is the one of the outgoing photons that remain 
for ever at the location of the light-like horizorP3. Combining Eqs. (|18fl and 
(|19fl we obtain m and S as functions of m lrT and uj : 



(20) 



y/l- 4mf„ OS 



2 t ,' 2 
irr 



S = , ■ (21) 

V 1 — 4m? r 

This is the right thing to do since lo is the b asic va riable describing each equi- 
librium configuration calculated numericall}EE21, and because the irreducible 
masses a re the o nes which are held constant along the numerical evolutionary 
sequenceJSEH. In the limit of small rotations we have 

S = I uj + O (u 3 ) , (22) 

where I = 4m 3 rl . is the moment of inertia of the Kerr black hold^l Next the 
total mass-energy is 

m = m irr + i/w 2 + 0(w 4 ) . (23) 

It involves the standard kinetic energy of the spin. 

To take into account all the spin effects our first task is to replace all the 
masses entering the energy function JSJ-© by their equivalent expressions 
in terms of uj and the two irreducible masses. It is clear that the leading 
contribution is that of the spin kinetic energy given by Eq. I|23|l and comes 
from the replacement of the rest mass-energy Mc 2 (where M = mi + m-i). 
From Eq. Q23J) this effect in the case of corotating binaries is of order uj 2 , which 
means by comparison with Eqs. that it is equivalent to an "orbital" 

effect at the 2PN order. Higher-order corrections in Eq. I|23|) . which behave 
at least like to 4 , will correspond at least to the orbital 5PN order and are 
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negligible for the present purpose. In addition there will be a subdominant 
contribution, of the order of w 8 / 3 equivalent to 3PN order, which comes from 
the replacement of the masses into the "Newtonian" part, proportional to 
x oc a; 2 / 3 , of the energy E [see Eq. (J2J)]. With the 3PN accuracy we do not 
need to replace the masses that enter into the post-Newtonian corrections in 
E, so these masses can be considered to be the irreducible ones in these terms. 

Our second task is to include the specific relativistic effects due to the 
spins, namely the spin-orbit (S.O.) interaction and the spin-spin (S.S.) one. In 
the case of spins Si and S2 aligned parallel to the orbital angular momentum 
(and right-handed with respect to the sense of motion) the S.O. energy reads 



Ei 



s.o. 



= -n (Muj) 5/3 



3M 2 



Si 



3 M 2 



So 



(24) 



Here we are employing the formula given by Kidder, Will and WisemarE^ 
who have computed the S.O. contribution as expressed by means of the or- 
bital frequency u> (this is what we need in view of the comparison with the 
numerical worl P ™ >). The basis of their computation is the work of Barker 
and O 'Conneli^l who obtained the formula g iven in terms of the orbital sep- 
aration r. The derivation of Eq. 1)24(1 in RefP^l takes correctly into account 
the fact that the relation between the orbital separation r (in a given coor- 
dinate system) and the frequency lu depends on the spins. We immediately 
infer from Eq. Il22[) that in the case of corotating black- holes the S.O. effect is 
equivalent to a 3PN orbital effect and thus must be retained with the present 
accuracy [with this approximation, the masses in Eq. I|24[l are the irreducible 
ones]. As for the S.S. interaction (still in the case of spins aligned with the 
orbital angular momentum) it is given by 



E s .s. =M ^(M W ) 2 ^% . (25) 

The S.S. effect can be neglected here because it is of the orbital order 5PN 
for corotating systems. 

Summaryzing, the contributions due to the corotating spins at the 3PN 
order are three in all : the main one is that of the spin kinetic energy and 
arises at the 2PN order; then we have two subdominant contributions at the 
3PN order coming respectively from a mixing between the spin kinetic energy 
and the Newtonian orbital energy, and from the S.O. interaction (|24p. We 
can neglect all the other terms. Summing up these three contributions we 
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find that the energy of the corotating spins is (coming back to the notation 
of Section 2) 



The complete 3PN energy of the corotating binary is the sum of Eqs. ©-© 
and (|26() . Notice that we must now understand all the masses in ©-(JSJl and 
(12611 as being the irreducible masses — we no longer indicate the superscripts 
"irr" — , which for the comparison with the numerical work must be assumed 
to stay constant when the binary evolves. 

The Figure ^ (issued from the wor presents our results for £"ico m 
the case of irrotational and corotational binaries. Since /\^E cmot , given by 
Eq. J2HJ), is at least of order 2PN, the result for lPN corot is the same as 
for 1PN in the irrotational case; then, obviously, 2PN corot takes into account 
only the leading 2PN corotation effect [i.e. the spin kinetic energy given by 
Eq. I|23(l ]. while 3PN corot involves also, in particular, the corotational S.O. 
coupling at the 3PN order. In addition we present in Figure ^ the numerical 
point obtained by numerical rel ativity under the assumptions of conformal 
flatness and of helical symmetrjESEll As we can see the 3PN points, and 
even the 2PN ones, are rather close to the numerical value. The fact that 
the 2PN and 3PN values are so close to each other is excellent, and confirms 
the good accuracy of the approximation we concluded in Section 3. In fact 
one might say that the role of the 3PN approximation is merely to "prove" 
the value already given by the 2PN one (but of course, had we not computed 
the 3PN term, we would not be able to trust very much the 2PN value). 
As expected, the best agreement we obtain is for the 3PN approximation 
and in the case of corotation : i.e. the point 3PN corot . However, the 1PN 
approximation is clearly not precise enough, but this is not surprising in the 
highly relativistic regime of the ICO. 

In conclusion, we find that the location of the ICO computed by numerical 
relativity, under the helical-symmetry and conformal-flatness approximations, 
is in good agreement with the post-Newtonian prediction. (See Ref.^2 f or the 
results calculated within the EOB method at the 3PN order, which are close to 
the ones reported in Figure^) This constitutes an appreciable improvement of 
the previous situation, because we recall that the ear lier es timates of the ICO 
in post-Newtonian theoryS2l and numerical relativitj@^^^ strongly disagreed 
with each other, and do not match with the present 3PN results (see Refill 
for further discussion). 

Notice in Figure ^ the difference between the energies of the corotational 
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(26) 



-0.01 



2pN corot 

% 3PN corot 
* NUMERICAL 
# 3PN 
* 2PN 



1 PN C ' 
1 PN 



0.3 

GM&j ir 



Figure 1. Results for the binging energy -Eico — Mc 2 versus uico i n the equal-mass case. 
The asterisk marks the result calculated by numerical relativity. The points indicated by 
1PN, 2PN and 3PN are computed from Eqs. and correspond to irrotational binaries. 

The points denoted by lPN corot , 2PN corot and 3PN corot come from the sum of Eqs. @J-0 
and 1261 . and describe corotational binaries. Both 3PN are 3PN corot are shown for u; s = 0. 



and irrotational configurations which amounts to approximately 3.5 x 10~ 3 
units of Mc 2 . We find that this amount is mainly made of the (positive) 
kinetic energy of the spins, which is about +5.0 x 10~ 3 , but that the S.O 
interaction at the 3PN order, equal to —1.3 x 10~ 3 , reduces slightly the effect, 
while the coupling between spin kinetic and orbital terms is quite negligible, 
of the order of — 0.2 x 10 ~ 3 . We feel that these numerical values are physically 
satisfying, and well in accordance with our order-of-magnitude estimates in 
Section 3 which showed that in the regime of the ICO the post-Newtonian 
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series "converges well", with the successive post-Newtonian approximations 
becoming numerically smaller and smaller. We have checked that the higher- 
order spin effects like the S.S. interaction [see Eq. Q25fl] which arises at the 
5PN order for corotating systems are completely negligible numerically (the 
S.S. effect is about +0.02 x 10" 3 ). 

Recently Damour et computed the effects of spin in corotating black 
hole binaries within the effective-one-body (EOB) approach. They find a 
result which differs from our computation presented above. Basically they 
obtain that the kinetic energy of the spins is about +5.0 x 10~ 3 , in agreement 
with ours, but that the S.O. effect is much larger than in our computation, and 
nearly comp ensa tes in their approach the spin kinetic energy. The net result 
given m RefEl for the total effect of spins is then about +1 x 10 3 , much 
smaller than our own result. The surprising fact (in the author's opinion) 
is that in the EOB method a relativistic effect like the S.O. coupling, which 
according to the a priori expectation based on the estimates of Section 3 
should be small, can nearly cancel out a "Newtonian" term, which one would 
physically consider to be the dominant one. Although Refs!""l an are 
in very good agreement for the ICO of corotational binaries, and in rather 
good agreement for irrotational binaries, their discrepancy concerning the 
numerical value of the S.O. effect remains unexplained. But let us try a 
guess. As we argued in Section 3 the resummation techniques are probably 
justified only when the post-Newtonian coefficients in the energy function get 
bigger and bigger as the order of approximation increases. Our guess would 
be that perhaps the EOB method is somehow forced to attribute in front 
of the S.O. term a big post-Newtonian coefficient in order to maintain its 
internal consistency. With an abnormally big coefficient one might explain 
the discrepancy. 
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